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Abstract 

We present the analytical calculation of the four-loop quark mass anomalous dimension 
in Quantum Chromodynamics within the minimal subtraction scheme. On the basis of 
this result we find that the so-called invariant quark mass is a very good reference mass for 
the accurate evolution of the running MS quark mass in phcnomcnological applications. 
We also obtain for the first time a complete 4-th order perturbative QCD expression for a 
physical quantity, the total Higgs boson decay rate into hadrons, and analyze the infrared 
fixed point for this case. 
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The renormalization group equations play an important role in our understanding of 
Quantum Chromodynamics and the strong interactions. The beta function and the quarks 
mass anomalous dimension are among the most prominent objects that appear in these renor- 
malization group equations. 

The quark mass anomalous dimension in QCD was calculated at the 3- loop level in [jy]. In 
a previous paper we calculated analytically the 4-loop QCD /3-function. Here we complete 
this 4-loop project by obtaining the quark mass anomalous dimension in the 4-loop order. 

We will apply the obtained anomalous dimension to show the relevance of the invariant 
quark mass for an accurate evolution of the running quark masses in phenomenological ap- 
plications. We also apply the 4-loop result to analyze a possible QCD infrared fixed point in 
the case of Higgs boson decay into hadrons. 

Throughout the calculations in this article we use the technique of Dimensional Regu- 
larization || for the regularization of ultraviolet divergences and the Minimal Subtraction 
(MS) scheme [|| or its standard modification, the MS scheme ||, for renormalization. The 
dimension of space-time is defined as D = 4 — 2e, where e is the regularzation parameter 
fixing the deviation of the space-time dimension from its physical value 4. 

Let us write the full Lagrangian of Quantum Chromodynamics with massive quarks in 
the covariant gauge 

1 n/ _ 

L = --G^ u G aflu + ^2ip q {Hp - m q ) ip q + Lgf + L gc 

q=l 

G 'fif = 9^A U — d v A^ + gf A^A U 
[D^ = Sydp-igAm (1) 

The gauge-fixing and gauge-compensating parts of the Lagrangian in the covariant gauge are 

L gc = d^ a {d^ a -gr bc JA^) (2) 

The fermion fields (the quark fields in QCD) tp q have a mass m q and transform as the fun- 
damental representation of a compact semi-simple Lie group, 5 = 1, ...,rif is the flavor index. 
The Yang-Mills fields (gluons in QCD) A^ transform as the adjoint representation of this 
group. uj a are the ghost fields, and £ is the gauge parameter of the covariant gauge. 

T a are the generators of the fundamental representation and f abc are the structure con- 
stants of the Lie algebra, 

rparpb rpbrpa ^ jabcrpc 

In the case of QCD we have the gauge group SU(3) but we will perform the calculation for 
an arbitrary compact semi-simple Lie group G. 

The beta function and the quark mass anomalous dimension are defined as 

am [i L 

= -(3 Q a 2 - f3 ia 3 - (3 2 a 4 - (3 3 a 5 + 0(a 6 ) (4) 
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din m q 
din (j? 



= ~7o« - 7i« 2 - 720 3 - 73a 4 + 0(a 5 ) 



(5) 



where m q = m q (fj?) is the renormalized (running) quark mass and /J, is the renormalization 
point in the MS scheme, a = a s /ir = g 2 /4-K 2 where g = g(fj> 2 ) is the renormalized strong 
coupling constant of the standard QCD Lagrangian of Eq. (|l]). 

This normalization for the basic perturbative expansion parameter a differs from the 
normalization g 2 /lQn 2 which naturally appears in multiloop calculations and was therefore 
used in Here we adopt the different normalization a = a s /ir because it is the commonly 
accepted normalization in applications. 

The 4-loop QCD /^-function in the MS scheme reads 



fa 
fa 
fa 
fa 

Or in a numerical form 

fa * 

fa « 
fa « 
fa * 



1 
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1 

64 
1 
256 



11 
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2857 



38 



-n 
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325 



2 18 
/ 149753 
V 6 



7 54 
+ 3564C 3 ) 



77 2 
"7 



/ 1078361 6508 
V 162 h 27 



/ 50065 6472 A 
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81 
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(6) 



2.75 - 0.166667n/ 
6.375 - 0.79166771/ 
22.3203 - 4.36892n/- + 0.0940394n 2 f 



'f 



114.23 - 27.1339™/ + 1.58238nJ + 0.0058567n 



(7) 



To calculate the quark mass anomalous dimension j m we need to calculate the renormal- 
ization constant Z m of the quark mass 

m B = Z m m 

where mg is the bare (unrenormalized) quark mass. Since the anomalous dimension 7 m does 
not depend on masses within the MS scheme it is formally the same for all quarks. We will 
therefore omit the flavor index q in m q when it is irrelevant. The expression of the 7 m -function 
via Z m is given by the following chain of equations 



d In tub 
din fi 2 







d In Z„ 



+ 



din m 



din fi 2 din fi 
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d\nZ m dlnZ m da dlnZ m fl , dZ$ 

* lm = W = ~^a~d^ = —da~ [ - £a + m = a ^a~ (8) 
where one uses the fact that the bare mass tub is invariant under the renormalization group 
transformations. Here da/ (din n 2 ) = —ea + (3(a) is the /3-function in L>-dimensions, Z^ is 
the coefficient of the first e-pole in Z m defined below. 

Renormalization constants within the MS-scheme do not depend on dimensional param- 
eters (masses, momenta) Q and have the following structure: 

i °° 7 {n) ( s 

n=l 

Since Z m does not depend explicitly on /i and m, the 7 m -function is the same in all MS-like 
schemes, i.e. within the class of renormalization schemes which differ by the shift of the 
parameter fi. That is why the 7 m -function is the same in the MS-scheme 0] and in the 
MS-scheme §]. 

We obtained the mass renormalization constant Z m by calculating the following two renor- 
malization constants of the Lagrangian: Z^ of the bilocal operator i^ip and Z^ of the quark 
field (i.e. of the inverted quark propagator). The connection between renormalized and bare 
quantities is defined as follows 

Z 



= Z l J 2 ^ (10) 

Then Z m = Z-x,/Z^. In order to calculate the 4- loop approximation of Z m we needed 
to calculate the 4-loop ultraviolet counterterms of the quark propagator and of the Green 
function 

= / dxdye^y(0\T{^(x)[^(y)i;(y)}m}\0) (11) 

These 4-loop ultraviolet counterterms were calculated using a technique which is based on 
the method of infrared rearrangement and is described in Ref. || . This general technique 
reduces the calculation of the counterterms to the direct calculation of 4-loop massive vacuum 
(bubble) integrals and provides a procedure that is well suited for the automatic evaluation of 
large numbers of Feynman diagrams. The calculations were done with the symbolic manipu- 
lation program FORM || and a specially designed database program that stores the results 
of the individual diagrams and adds them in the end. For the present calculation we needed 
to evaluate of the order of 10000 4-loop diagrams. These diagrams were generated with the 
program QGRAF §. 

We obtained in this way the following result for the 4-loop 7 m -function in the MS-scheme 

70 = \[3C F ] 

71 = — \-Cp + —CfCa — —CpTprif 

16 2 F 6 3 ' 
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(12) 



Here C is the Riemann zeta function (C3 = 1.2020569 •• •, C4 = 1.0823232- •• and C5 = 
1.0369277- • •). [T a T a ]ij = C F 5ij and f acd f bcd = C A 5 ab are the quadratic Casimir operators 
of the fundamental and the adjoint representation of the Lie algebra and tv(T a T b ) = T F 5 ab 
is the trace normalization of the fundamental representation. A^^ is the dimension of the 
fermion representation (i.e. the number of quark colours) and nj is the number of quark 
flavors. We expressed the higher order group invariants in terms of contractions between the 
following fully symmetrical tensors: 

1. 



d abcd 



mbcd 

d A 
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-Tr 



rparpbrpCrpCl _j_ rparpOrpdrpC _|_ rparpCrpOrp, 



drpC 



-<arpCrpbrpd 



rparj-iCrj-idrpb _|_ rparpdrpbrj-iC _|_ rparpdrpCrpb 
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(J a (J^ (J C (jd _|_ (J 0, (J^ (J^ (J C -|- (J* 1 (J C (J^ (J° 



+C a c c c d c b + C a c d c b c c + C a c d c c c b 



(13) 



(14) 



where the matrices [C a ]fc c = —if are the generators in the adjoint representation. 



The result of Eq. (O) is valid for an arbitrary semi-simple compact Lie group. The result 



for QED (i.e. the group U(l)) is included in Eq. (|T|) by substituting C A = 0, d a A bcd = 0, 
C F = 1, T F = 1, (df cd ) 2 = 1, Np = 1. The n 3 f and nj terms in this 4doop result for QED 
agree with the literature |j~C| ] where the leading and next-todeading large-n/ terms for the 
QED gamma-function were calculated in all orders of the coupling constant. Furthermore 
it is interesting to note that for the choice C F = C A = T F ,nf = 1/2 and dp = d A , which 
corresponds to the case of N=l super symmetry, all ( terms and terms with the tensors d 
cancel. This is analogous to the case of the 4-loop /3-function M. 
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The result of Eq.([T^) was obtained in an arbitrary covariant gauge for the gluon field. 
This means that we have kept the gauge parameter £ that appears in the gluon propagator 
i[—gt Ml/ + (1 — Qq^q" / (q 2 + ie)]/(q 2 + ie) as a free parameter in the calculations. The explicit 
cancellation of the gauge dependence in 7 m gives an important check of the results. The 
results for individual diagrams that contribute to 7 m also contain several constants specific 
for massive vacuum integrals. The cancellation of these constants at various stages in the 
calculation provides additional checks of the result. We note that at the 4-loop level £4 and 
£5 (Riemann zeta function of arguments 4 and 5) appear as new constants. 

For the standard normalization of the SU(iV) generators we find the following expressions 
for the color factors (more details about the color factors for various other groups can be 
found in 01) 



1 



Tp = — 
2 



jabcd^abcd 



N F = N, 

(N 2 -l)(N 2 + 6) 
48 ' 



C A = N, C F 



N 2 



1 



jabcd^abcd 
F F 

N F 



2N ' 

(iV 2 -l)(iV 4 -6iV 2 + 18) 
96iV 3 



Substitution of these color factors for N = 3 into Eq. (12) yields the following result for QCD 
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/ 



(16) 



Let us now consider the solution of the evolution (i.e. renormalization group) equation 
Eq.(|5|) for the quark mass 



m g (/i 2 ) = m g (/io)exp 



( a ^ 2) da nm(a') 

a(^) a Pip!) 



m q exp 



" a m 



(17) 



where we formally define the renormalization group invariant (i.e. independent of /j, 2 ) quark 



mass m„ as 



fh q = m q (i4)exp ^- J 



da WY 



(18) 
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One can expand Eq.(|l7|) in a to obtain the following perturbative solution to the evolution 
equation 



m q {ii 2 ) 



7 7o/Aj 



1 + A x a + (A\ + A 2 ) y + ( 
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+ 0(a 4 ) 
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where 
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This gives us the following 4- loop expansions for the running MS quark masses m q (p z 
for quark flavors q=s,c,b,t 
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7T 
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7T 



(20) 



where a s = a s (/i 2 /A|^). 

In phenomenological applications one uses the MS quark masses at a characteristic scale 
H of a considered process to sum large logarithms. From the expansions (|20| ) we conclude 
that the invariant mass rh q is a good reference mass for the accurate evolution of the MS 
quark masses to the necessary scale /i in phenomenological applications. The perturbative 



coefficients of the solutions (20) of evolution equations for the MS quark masses are small up to 
and including the 4-loop level. On the contrary, it is known that the equations expressing the 
MS quark masses m q via the pole quark masses M q have worse convergence of the perturbative 



expansions [11]; for a review see Ref. [12]. We should emphasize in this respect that the 
invariant mass is a fundamental gauge-invariant object that naturally appears in the solution 
of the renormalization group equation. 
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The analysis of the QCD infrared fixed point 



Let us apply the obtained 4-loop quark mass anomalous dimension to the analysis of the QCD 
infrared fixed point in the case of the Higgs boson decay into hadrons. The infrared fixed 
point was analyzed at the third-order of perturbative QCD in Ref. K3| for the total hadronic 



Higgs decay and in [14, 15] for electron-positron annihilation into hadrons. Recently the total 
hadronic decay width of the Higgs boson was calculated at the 4-loop level of perturbative 
QCD |1(| in the limit of massless quarks. 
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1 



vr - -Cs 

162 3 
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4V2tt 

+ (164.1392 - 25.7711971/ + 0.258974n}) a 

_ a s (M H ) 



(21) 



m q = m q (MH) and Mjj is the Higgs mass. We should note that if the 
Higgs boson is lighter than the top quark, then additional contributions appear from the 
so-called singlet diagrams due to non-decoupling of the heavy top quark in this channel [17], 
see also Ref. [18]. We will neglect these contributions. 

Together with the 4-loop quark mass anomalous dimension it allows for the first time to 
perform the analysis of the physical quantity at the fourth-order of perturbative QCD, i.e. 
with four known perturbative QCD a s -terms. 

For the analysis of the infrared fixed point it is convenient instead of Th to introduce the 
function 

_ldhx{V H /M H ) _ ^ l ar y d]xiT(a) 
~2~ 
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And the full result for R becomes 
11 

n = 12 - — n } 

« 12-0.6111171/ (24) 

7189 209 , 429, / 2995 5 2 17 

ro = vr G + ^f 1 Ti" H Cs 

2 36 48 8 s ' V 144 12 4 w 
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376.12 - 135.72n/ + 7.2045n}-0.05895n) (26) 

One may notice that the C4-te-rms coming separately from the 4-loop mass anomalous dimen- 
sion and the 4-loop decay width cancel in the final expression for R(a). This cancellation of 
^4-terms provides an extra cross-check for both the calculations of 7 m and Th- The remain- 
ing 7r 2 -terms come from the imaginary parts of logarithms in the process of the analytical 
continuation from the Euclidean to the Minkowski region. 

For the analysis of the infrared fixed point we will apply the approach of effective charges 



19]. This approach to the scheme dependence of perturbative series is known, see e.g. 
Ref. |14| , to give values for physical quantities close to another distinguished renormaliza- 
tion scheme - optimized perturbation theory p0[. Within the approach of effective charges 



one defines the whole calculated perturbative expansion as a new effective charge, i.e. 

R(a) = a + via 2 + r 2 a 3 + r 3 a 4 = a c g (27) 

The evolution equation for the effective charge (or equivalently for the physical quantity itself 
in the given order of perturbation theory) is governed by an effective /3-function 



dR 



d\nM 2 H 



MR) 

-P0R 2 - P1R 3 - PfR 4 - PfR 5 + 0(R 6 ) (28) 



Thus in the effective charge approach one deals with only one asymptotic perturbative series 
for f3 e ff(R) instead of two series for R(a) and for (3(a) in the MS-scheme. 

The third and the fourth coefficients of the effective (3- function (which are scheme-invariant 
since they govern the evolution of the physical quantity R(a)) are expressed via the MS 
coefficients as follows 



9 



pf = - Pm + Po(r 2 - rf) 

Pf = p 3 -2p 2 r 1 +p 1 r 2 1 +2p (r 3 -3r 1 r 2 + 2rf) (29) 
The numerical results are 

Pf ps -196.464 + 26.1453^- 1.38317^ + 0.0290035^ (30) 

pf ps 3305.698 - 700.571ra/ + 65.1914r^- 2.89465n^ + 0.0499219n) (31) 

In the third order of QCD this /3 e fj-function has a positive zero R.W « 0.15 Hi which 



practically does not depend on the number of quark flavors for nj = 3,4,5,6. In Ref. [13] it 
was considered to be a spurious fixed point, but it could have been interpreted as a real QCD 
infrared fixed point which would indicate the possibility of applying perturbation theory till 
zero energy. This is why it is important to check the stability of this zero under the 
inclusion of the fourth order of perturbative QCD. One can see from Eq. ( pi] ) that at the 
fourth order level the positive zero of the effective beta function disappears. This supports 
the idea that the infrared fixed point for the Higgs boson decay into hadrons at the third 
order of perturbative QCD is indeed a spurious fixed point. 
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Just before completing this paper a similar paper appeared [21] with the mass anomalous 



dimension to four loops. It agrees completely with our formulae |15| and 
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